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Abstract 

Disoriented chiral condensates or long wavelength pionic oscillations and 
their interaction with the thermal environment can be a significant source of 
dileptons. We calculate the yield of such dilepton production within the linear 
sigma model, both in a quantal mean-field treatment and in a semi-classical 
approximation. We then illustrate the basic features of the dilepton spectrum 
in a schematic model. We find that dilepton yield with invariant mass near 
and below 2m,r due to the soft pion modes can be up to two orders of magni- 
tude larger than the corresponding equilibrium yield. 
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1 Introduction 



If the quark masses are neglected in two-flavor QCD, there would be a continuum of 
degenerate vacua in which chiral symmetry is spontaneously broken. These vacua are 
mutually related by SU{2) rotations and can be characterized by non-vanishing scalar 
(sigma), {a) = (qq), and pseudo-scalar (pion), (tTj) = {q'j5Tiq), quark condensates. 

In the "Baked Alaska" scenario proposed by Bjorken, Kowalski, and Taylor a space- 
time region protected by a hot shell from the normal vacuum (where (tTj) = 0, (a) ^ 0) 
can relax to a misaligned vacuum with a non-vanishing pion condensate, or disoriented 
chiral condensate (DCC). After the dispersion of the hot shell, the DCC will then decay 
into normal vacuum by coherent emission of soft pions. In relativistic heavy-ion collisions, 
if the system maintains quasi-equilibrium through the chiral phase transition, the quark 
masses, though small, will prevent the correlation length from growing indefinitely. How- 
ever, recent numerical studies of the linear sigma model 0, |], ||, |, ||, |^, ^, |[ have shown 
that a rapid cooling-like quenching can drive the system far out of equilibrium and lead 
to significant amplification of the soft pionic modes. The resulting occupation numbers 
may then become large and lead to the emission of many pions in the same isospin state. 
In such an ideal scenario, the neutral pion fraction / exhibits an anomalous distribution, 
P{f) = 1/2^7) which has been suggested as an experimental signal [|10], [11], |12|. However, 



if several separate domains are formed {i.e. if the size of the system is large in comparison 
with the correlation length), as may well occur in heavy- ion collisions, the signal is corre- 
spondingly degraded and the distribution approaches its normal form, an approximately 
normal distribution centered around /=| 0, [Hj- More advanced methods of analysis 
would then be needed, such as a the use of wavelets or cumulative moments . 
In addition to the hadronic signals, the electro-magnetic DCC signatures have also 



been also addressed [|T7|, |T8[. Since the electro- magnetic current is given by the third 
component of the isovector current, the isospin oscillation of a coherent pion field may 
produce photons and dileptons which could provide information on the early dynamical 
evolution of the DCC. In Ref. Jl^, this emission from the oscillation of the coherent pion 
field was found to be sizable only for very small invariant masses, and thus it is hard to 
observe experimentally due to the large background from ttq Dalitz decays. However, a 
DCC may also contribute to the incoherent production of dileptons. In particular, high- 
momentum incoherent pions may annihilate with the coherent DCC pions. Since the pion 
phase-space density in a DCC is comparatively large and well localized in momentum 
space, this process should lead to a considerable enhancement in the dilepton spectrum 
at finite invariant masses (~ 2171^^) that is rather narrow in invariant mass as well as 
transverse momentum. 

Since the dynamical evolution leading to DCC formation is far from equilibrium, 
we shall study the electromagnetic production processes within the framework of time- 
dependent field theory, using both a quantal mean-field treatment |jl9| and a semi-classical 
approximation As these calculations treat coherent and incoherent production on an 
equal footing, the results do not depend on how that distinction is made. Therefore, to 
provide a simple understanding of the phenomenon, including the different contributions 
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from coherent and incoherent production, we shall subsequently illustrate the essential 
features by means of a schematic model. 

This paper is then organized as follows. In the next section, we first briefly review the 
formulas needed to calculate the dilepton production. We then describe the dynamical 
calculations obtained with the linear sigma model and present the corresponding results 
for the dilepton yields. Finally, we turn to the schematic model. 



2 Dilepton production from DCC 

In general, in a non-equilibrium system, such as the dynamical evolution of DCC fields, 
the in and out states are not asymptotic states and the density matrices describing the 
system are not diagonal (except in isospin and charge). In this paper, we will neglect 
quantum effect caused by the off-diagonal matrices of the system in the calculation of 
the dilepton production. Therefore, the dilepton production yield can be given by the S 
matrices of the electromagnetic transition between different states , 



lA id —in-i T — ii\n-\ r ^ ^- / 



where we have summed over the final states, Z = Tr[p], and 

5 =(1-^)^/^(1 + ^). (2) 

Since we are interested in the production of electron-positron pairs, we shall neglect the 
lepton mass mi in the following, i.e. B = 1. 

We have modified the temporal integration boundaries from the asymptotic times 
t = ±cx) to finite initial and final times ti and tf. This is suitable for an initial- value 
problem where the initial conditions are fixed in non-asymptotic states which are assumed 
to be formed in a relativistic nuclear collision. 

The Lagrangian density of the linear sigma model is given by 

C = \d,<pd^<P - ^A(02 - v')' + Ha , (3) 

where (p = (cr, tt) are the chiral fields in 0(4) representation. The parameters. A, v and 
H are determined by the pion and sigma mass, 171^^, m^, and the pion decay constant /^r. 
In order to include electromagnetic (EM) processes, we can introduce an EM field in the 
charged sector of the above Lagrangian. The symmetrized EM current density coincides 
with the third component of the isovector current density and is given by 

j/.(a;) = ^[7r"^(a;) d,, T^ix) - n{x) ti\x)] = Tii{x)d^TX2{y) - 7r2{x)d^7ii{y) , (4) 
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where the complex charged pion field operators are related to the Cartesian components 

by ^ ^ 

7r(x) = --^[ni{x) + i7i2ix)] , n\x) = --^[7ri(x) - in2{x)] . (5) 

We shall neglect the quadratic coupling in the gauged linear sigma model and the anoma- 
lous electromagnetic coupling of 7r°, which all contribute to the dilepton production only 
to higher orders in the fine structure constant a = e^/An. 



2.1 Mean- field treatment 

In the mean-field treatment, the four-point functions in the current-current correlator in 
Eq. (|l|) are given as products of two- and one-point functions similar to the decomposi- 
tion of a free field four-point function according to Wick's theorem. Since the one-point 
functions then represent the average field they vanish for the pions. For uniform density 
matrices the contributions from two-point functions in the coincidence limit vanish since 
= 0. Furthermore correlators of the type {Ti{x)Ti{y)) or (7r^(x)7r^(?/)) contribute 
only to zero-momentum processes and and vanish altogether for the initial conditions 
chosen in the present calculations. Therefore, the only remaining components for the 
current-current correlator are of the type {'n{x)7['^ {y)) . The current-current correlator 
then takes the form. 



W^^{x,y) = {TT^{x)n{y)){d^fi{x)d^7i^{y)) + {d^Tf^{x)d^Tf{y)){7r{x)7i^{y)) 

- {d^n{x)n\y)){n\x)d^7r{y)) - {df,7r\x)n{y)){n{x)d^n^{y)). (6) 

For a system in thermal equilibrium at the temperature T, the ensemble average is 
{■)th = Tr[/3t/i-], where pth ~ e~^^^. For free pions, the eigenstates are plane waves and 
the thermal ensemble average can be expressed accordingly, 

where nf = l/{e^^'^ — 1) are the Bose-Einstein occupation numbers for the charged pions. 
The hadronic tensor (|^) entering into the dilepton-production expression then becomes 

- ' :(27r) |(A;i - k2)^,{l^i - k2)„ 



d^x J 2cJi(27r)3 2LJ2(27r)3' 

■ [n~unl^\q 'h- ^2) + (1 + )(1 + 05^(g + h + /c^)] 
+ (A;i + fc2)M(^i + ^2).K(l + %J+<(l + 0]5'(g-fci + fc2)} , (8) 

where we have taken advantage of the translational invariance in space and time to obtain 
the ensemble average rate of production by dividing by the corresponding four-volume 
Vitj-U). 
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The three terms in the above relation correspond to pion pair annihilation, creation 
and bremsstrahlung, respectively. The factor 1 + is a result of the Bose enhancement 
in the final states. For dilepton production (g^ > 0, go > 0) only the pair annihilation 
term remains due to energy and momentum conservation. However, as we shall discuss 
briefly in section |^, once the pion dispersion relation develops several branches due to 
interactions, one can simply replace (1 + nf) by Y^i^k^i) + (^))] iii fhe above 

equation with summation over the number of branches in the dispersion relation. Since 
dileptons can be emitted via the transition from one branch to another, all the terms, 
except the pair creation process, are allowed. 

When the system is not in equilibrium, as is generally the case in a collision scenario, 
the current-current correlator cannot be expressed explicitly. However, in the present 
mean-field approximation the pion field vanishes on the average, and so the pion field 
operators can be expanded in terms of the annihilation and creation operators for medium- 
modified charged pions, 

^ J (0je^'-^[/,(t)a, + mhU] . (9) 

where the coefficients are the time-dependent mode functions fk(t). 

The time evolution of the mode functions can be obtained in the mean-field approxi- 
mation and we can then obtain the two-point correlation functions in Eq. (P), as described 
in Ref. J^, 19|. To solve the time evolution of the mode functions one needs to specify 



the initial conditions of the mean-field (cr) and its time derivative, the mode functions 
and their time derivatives, the average occupation numbers for each mode {a\ak) = n'l, 
(blbk) = and their pair correlation (cikb-k) = Fk- 

In order to have a finite set of renormalized equations, we have to choose the mode 
functions so that the high-momentum modes coincide with the zeroth order adiabatic 
vacuum described by 



^2uJk{to) 



-iujkito) - 



2uJkito) 



fkito) , (10) 



with uKto) = < x{to) > and X = A(0V8 - v^). 

In this study, we will select ensembles with n'^ = ra^ and vanishing pair correlations, 
Fk = 0. By calculating two- and four-point functions the parameters A, H, and v were 
chosen to give the best fit to the physical observables /^r, m^r, and the phase shifts 
1^. In these simulations the bare coupling constant is A = 20, the momentum cutoff is 
A = 1 GeV, the pion mass is 135 MeV and the pion decay constant is /^r = 92.5 MeV. The 
resulting mass of the a is 340 MeV. 

We investigate two different initial ensembles and evolve them from time ti = —50 
fm/c to t/ = 50 fm/c. All cases were prepared with equivalent energy densities and with 
zero kinetic energy density associated with the a mean-field. The vacuum adiabatic mode 
functions are assumed as 

fk{to) = , fk{to) = -i^kfkito) . (11) 
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The first one is an ensemble which at the initial time {t = —50 fm/c) has a mean 
a field which is slightly perturbed around a thermal ensemble prepared at T=100 MeV. 
The initial quasiparticle frequency in this case is given by 

ujI = e+<x> , (12) 

with the initial value < x >= (138MeV)^. The initial value of the mean a field is 
87 MeV. The average quasiparticle population is given by the corresponding Bose- 
Einstein distribution and the energy density is then e = 13.7MeV/fm^. At this relatively 
low temperature, the ensemble should be close to that of free pions and the results may 
thus serve to assess the quality of the calculation. Subsequently, we will call these initial 
conditions thermal. 

The second ensemble is prepared with a negative value of x order to emulate a 
quench scenarios. In order to initialize the unstable pion modes (those having a negative 
value of ujI), we have employed the relation ul = k'^ + + (< X > —f^i'^) exp {—f^/k'^), 
using the initial value < x >= —(132 Me V)^. [We note that at a large momentum the 
frequency agrees with the adiabatic vacuum frequency, as required by renormalization.] 
The initial value of the mean field is < cr >= 31 MeV and the energy density is e = 
13.6 MeV/fm'^, i.e. similar to the equilibrium scenario. In the quench case, the initial 
quasiparticle modes are unoccupied, = 0, corresponding to the local vacuum. 

Let us now turn to the results for dilepton production. In Fig. |1| we show the dilepton 
yield obtained by using the mode functions generated from the thermal initial conditions 
to calculate the current-current correlation function in Eq. (^ (full lines). Note that we 
have divided the yield by the total space-time volume, VAt, it is therefore the averaged 
production rate during the time interval At = tj — ti. We show the results for two different 
values of the three momentum of the dilepton, g = 50 MeV and q = 200 MeV. This is 
compared with the analytical result for a free-pion gas (dashed lines), 

1 dN^l\. _ go 4m^ f-^+ du ^ _ 
VTdMd^q AStc^M^ M^'L- g'^'^^"^' ^'^^ 

with cj^ = (go ± q^Jl — 4m^/M2)/2 and go = y/M^~+~q^. As compared to the analytical 
result for the free pion gas, the numerical calculation based on the thermal mode functions 
give a contribution below the two-pion threshold as well. This contribution is spurious as 
it arises from the finite resolution of the energy 6 function in our calculation: Since the 
Fourier transform in time extents only over a finite time interval, the energy is not perfectly 
conserved. As a result the terms proportional to ^^^(1 + n^^) Eq. (|^) do contribute to 
the time-like sector while with perfect energy conservation their contribution is restricted 
to space-like photons only. This spill-over from the space-like region is amplified because 
of the factor 1/g^ in Eq. (p. In addition to leading order this contribution depends on 
linearly on the distribution function Uk and is thus further enhanced by l/uk ^ 1 with 
respect to the pion annihilation term. 

In order to verify this, we have also performed a calculation with the mode functions 
of a free pion gas. The result is indistinguishable from the one obtained with the dynam- 
ical mode functions, showing that the interaction plays a negligible role at T=100 MeV. 
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By inspecting the individual contributions, we have also verified that the spurious con- 
tribution is indeed due to the terms proportional to nk^{l + rik^)- At the present stage, 
this spurious contribution forms an numerical background. It is also responsible for the 
slight disagreement between the analytical and numerical results for the free pions above 
the two-pion threshold. With this numerical background in mind, let us turn now to the 




Figure 1: Dilepton invariant mass spectrum for the thermal initial conditions for two 
values of the three momentum q (full lines). Also shown are the theoretical results for a 
free pion gas (dashed lines). 

comparison between the thermal and quench initial conditions, as illustrated in Fig. |^. 
Clearly, around an invariant mass of M ~ 2m^ the dilepton production rate from the 
quench initial conditions is about two orders of magnitude larger than the thermal pro- 
duction. This enhancement is less for larger momenta, refiecting the narrow momentum 
distribution of pions in the DCC fields. Furthermore, at low invariant masses there seems 
to be an enhancement as predicted in Ref. fl^. Because of our numerical background. 



we are not able to give a quantitative estimate of this enhancement. However, due to 
the strong background from the ttq Dalitz decay, this enhancement is difficult to measure 
experimentally. Therefore, let us concentrate on the enhancement around M ~ 2771^^. As 
expected from the narrow momentum distribution of the DCC pions, this enhancement 
is localized in invariant mass as well as in transverse momentum. The latter can be seen 
from Fig. ^ where the dilepton yield is plotted as a function of the three-momentum of 
the dilepton for an invariant mass of M = 300 MeV. The enhancement is confined to low 
momenta, below q ^ 300 MeV. While the localization of the enhancement in invariant 



6 



10-^ 




-IQ-12 I , , . , . , , , . , , , , , . 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
Dilepton mass (MeV) 



Figure 2: Dilepton invariant mass spectra for thermal (full lines) and quench (dashed lines) 
initial conditions. Shown are the spectra for two different values of the three dilepton three 
momentum q. 



mass nicely reflects the enhancement of low-momentum pion modes, it may be very diffi- 
cult to observe experimentally. At low transverse momentum the background from false 
pairs is largest and it remains to be seen if present detector designs allow a sufficiently 
accurate subtraction to make the extraction of this signal feasible. In order to get an 
idea about the experimental constraints, we have subjected our results to the CERES 



acceptance cuts |2^. The resulting invariant mass spectrum is virtually indistinguishable 
from the thermal one. Only if the momentum cuts can be pushed down to ^ 100 MeV 
will an enhancement of about a factor of ten remain. 



2.2 Semi-classical treatment 

In order to verify that our results are robust with respect to model details, we have 
also calculated dilepton production with a different treatment of the linear sigma model. 



namely the semi-classical approach described in Ref. pO|. In that treatment, the system 
is described by the real field 0(r, t) = (cr, tt), which is evolved by the classical equation of 
motion, 

[a + X{(f)^ -v^)]cj) = Ha , (14) 

where 0^ = o and a denotes a unit vector along the a axis. In the present context, 
when we are concerned with dilepton production in a macroscopically uniform system. 



7 




Dilepton momentum q (MeV) 

Figure 3: Momentum spectra for thermal (full lines) and quench (dashed lines) initial 
conditions for dilepton pairs of invariant mass M = 300 MeV. 

it is natural to enclose the system in a (sufficiently large) box with periodic boundary 
conditions. The spatial average of the field can then be considered as the order parameter, 
(^(t)=<(^(r, t)>, and the residual fiuctuations represent quasi-particle excitations relative 
to that constant field, 5(/)(r, t)=(/)(r, t) — (p{t). A similar separation holds for the time 
derivative of the field. 

In order to solve the above equation of motion, it is necessary to specify the initial 
field configuration, 0(r, tj), and the associated time derivative. Employing the method 
developed in Ref. [^, we sample these initial conditions from a thermal ensemble at a 
specified temperature T. The occupation numbers rik of the quasi-particle modes are 
drawn from the appropriate Bose-Einstein distributions which ensures that the initial 
state displays the proper quantum-statistical features. The resulting treatment is then 
akin to the Vlasov model often employed in nuclear dynamics at moderate energies, which 
is the semi-classical analogue of the Time-Dependent Hartree-Fock description. 

In order to make contact with the quantum mean-field studies described earlier, we 
consider two scenarios having the same energy, namely thermal equilibrium at T= 140 MeV 
and a corresponding quenched initial condition. The quenched ensemble has been gener- 
ated from the thermal ensemble by suppressing the field fiuctuations by a factor of 100 
and resetting the order parameter to 4'{ti) = (cro,0), with a vanishing time derivative, 
i/'(r, tj) = (0,0). Using o"o = 32 MeV ensures that the energy of the quenched scenario 
matches the thermal value, which is about 40 MeV/fm^. 

The numerical solution of the equations of motion of the chiral fields yields the evolu- 
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tion of the Cartesian components of the pion field, 7r(r, t), in addition to the sigma field 
a{r,t). The corresponding electromagnetic current density is then easy to extract, 



J^{x) = 7Tl{x)d^TT2{x) - 7l2{x)d^TCi{x) . (15) 

In order to calculate the associated invariant differential dilepton yield, we employ the 
following expression |pT| , 
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d^y J^ix) e-"i^--y^ My) (16) 

(17) 

This expression ignores the final-state Bose enhancement factors l+n^ which is justified 
when the occupation number n-^ is small, as is typically the case in equilibrium ||20|. In 



particular, for T=140 MeV the occupancy of the lowest quasi-particle mode is about 0.30. 
However, the occupancy may reach several units in the corresponding quench scenario and 
the formula (^) may then significantly underestimate the contribution from the softest 
pion modesQ. 

The last expression (|l^) recasts the dilepton yield in terms of the four- dimensional 
Fourier transform of the charge current, 

J^(g) = J d'^x J^{x) e^«^ = e*'^"*-^'^ '" . (18) 

The factorized form (|1^ is of great numerical convenience, since the transform J^{q) can 
be readily accumulated in the course of each dynamical history, and the contraction in 
(p!7|) need only be carried out at the end of the evolution. In (|1^) and (|I8]) the time 
integration extends over the duration of the observation, from tj to tf. In the equilibrium 
scenario, the yield then becomes proportional to At = — as well as to the volume V , 
and so a division by the four-volume VAt yields the corresponding invariant production 
rate d^N / {d^qiV /\t)). 

In general, we consider an entire sample of A/" individual evolutions, (r, t)}, where 
the label n enumerates the individual "events" in the sample. The resulting ensemble- 
average dilepton yield is then 

d^N 1 ^ rf^^W 

where d^N^'^'^/d'^ q is the contribution from the particular event n, obtained as described 
above. Since we consider ensembles that have translational symmetry, the current-current 

^ It may be noted that if the calculated pion field 7r(r, t) is assumed to represent a standard coherent 
state, then the quantal evaluation of the dilepton radiation rate would lead to the above expression ( |l6|) 
when the commutator terms are ignored; if those commutator terms were retained, then the final-state 
Bose enhancement factors would be recovered. 
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correlation function, -< J^(x) Jij{y) >~, will depend only on the spatial separation. More- 
over, in the special case of an equilibrium ensemble, its temporal dependence in equi- 
librium is only via the time difference. [In practice, the considered dilepton observables 
vary little from event to event, because the system is larger than the correlation length, 
and therefore sufficiently accurate results can be obtained on the basis of rather small 
samples.] 

In order to verify that the adopted method indeed leads to physically reasonable 
results, let us consider the production of back-to-back dileptons from a thermal gas of free 
pions. In that special case, the four-momentum of the dilepton is of the form q = (M, 0) 
and we are interested in masses M above 2171^^- The current-current contraction in (|l6|) is 
then especially simple and it is elementary to show that its ensemble average is given by 

-< M^) - 9"'^ My) ^ = ^ Ji^) ■ Jiy) >- = 2|vcp - 2CAc , (20) 

where we have employed the thermal correlation function of the charged pion fields, 

7ii{x)7ri{y) y = -< H2{x)7i2{y) y = C(r,t) = — ^ — cos(k ■ r - , (21) 

with being the thermal occupancy, hk = l/{exp{u!k/T) — 1), and (r,t) denoting the 
difference x-y. [The adopted sampling procedure ensures that the numerically extracted 
correlation function indeed yields this expression Since the back-to-back dileptons 



have vanishing momentum, q = 0, the Fourier transform over the separation r reduces to 
a spatial average and we readily find 

/ dv ^J{x)-Jiy)y= 4/ rfr|VC(r,t)|2 = 1 y [1 + cos2cUfet] . (22) 

Jv Jv k 

The remaining Fourier transformation over the temporal difference then restricts the 
contributions in the sum to those modes that have frequencies Uk near half the dilepton 
mass, M/2. Thus, in the continuum limit, when both the box and the time interval 
are large, we recover exactly the usual expression for production of dileptons by pion 
annihilation, 



^ - -(^? f^fdt [1 + cos2c.,tl e^^* - (i_^y ^23) 

d^qd^x ~ 3n^2J J (27r)3 J ^* u;l + ^ " 3 {2n)^ ) '^^^^ 

where no denotes the occupancy of pion states with the matching frequency ujq = M/2. 
Thus, at the formal level, the semi-classical method is indeed physically reasonable. 

In order to illustrate the numerical reliability of the calculations, we consider again the 
simple case of a free pion gas in thermal equilibrium, in analogy with what was shown in 
Fig. m. Figure ^ shows both the analytical result for the invariant dilepton production rate 
and the corresponding results obtained by integrating the equation of motion for the free 
fields, with the initial field configurations having been sampled from the associated thermal 
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ensemble as described above. As in the case of the mean-field treatment, the numerical 
calculation yields a reasonable result for dilepton masses above the pion annihilation 
threshold, but exhibit a divergent behavior for lower masses as a refiection of the pole at 
q = 0. The fiuctuations in the numerical results are a consequence of the much coarser 
grid employed in the semi-classical calculations: the system is confined within a cubic 
torus of size L=24 fm (with a grid spacing of 0.24 fm) and 41 modes have been included 
in each of the three Cartesian directions, for a total of 68,921 modes. Moreover, the time 
evolution has been performed only up to 50 fm/c. 




Figure 4: The invariant dilepton production rate in a gas of free pions in thermal equilib- 
rium at T = 100 MeV. The solid curves are obtained by solving the free-field equations 
numerically, while the dashed curves represent the corresponding analytical results. The 
heavy curves are for back-to-back production, q = 0, while the light curves are for a finite 
momentum of the dilepton, q = 300 MeV. 

We now move on to discuss the results of the numerical simulations in the two scenarios 
described above. Figure ^ shows the invariant production rate as a function of the dilepton 
mass M = -y/g^, for various magnitudes of its momentum q, in a display similar to Fig. 
The results are qualitatively similar to those obtained with the mean-field treatment: the 
quench scenario leads to a large enhancement around 400 MeV. In the present case the 
enhancement is only about one order of magnitude, because the coarser grid considered 
leads to a quicker damping of the oscillations in the order parameter and, consequently, 
the amplification of the soft pions modes is less extreme. Still, there is clearly a very 
significant effect of the non-equilibrium evolution. 
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Figure 5: The invariant mass spectrum for both thermal initial conditions (solid curves) 
and the corresponding quench scenario (dashed curves), for either back-to-back dileptons 
having q=0 (heavy curves) and dileptons with a finite momentum of q=150 Me V (light 
curves) . 

Finally, Fig. ^ shows the effective production rate d^N/{dMd^q (VAt)) as a function 
of the magnitude of the dilepton momentum q, for dilepton masses near M=300 MeV. 
Again, we see how the non-equilibrium evolution following the quench caused a large 
enhancement of the slow-moving dileptons. The results of the two calculations agree 
qualitatively: The quench initial conditions lead to a substantial enhancement of the 
dilepton yield at invariant masses close to M ~ 2771^^. The calculation based on the quan- 
tal mean-field approximation (section |2.1|) predicts a considerably stronger enhancement 
which is also more narrow in invariant mass as well as momentum. As we will explain 
in the following section this implies that in the quantal mean-field approximation the 
enhancement of the low momentum modes is stronger and narrower in momentum space 
than in the semi-classical approximation. Translated into the DCC language, the DCC 
fields generated in the quantal mean-field approximation have a larger field strength and 
are of larger spatial extent. There are several possible reasons for this quantitative differ- 
ence. First, the Bose enhancement factors, present in the mean-field treatment, certainly 
will increase the dilepton yield to some extent. Second, the comparatively coarse grid em- 
ployed in the semi-classical treatment leads to a faster dampening of the oscillations in the 
order parameter, resulting in a smaller amplification of the soft pion modes. Finally, the 
semi-classical treatment incorporates the mode mixing arising from the non-linear form 
of the interaction. This mechanism helps to equilibrate the system and thus redistribute 
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Figure 6: The production rate {d'^N/dMd^q {VAt)) as a function of the magnitude of 
the dilepton momentum q, for dilepton masses near M=300 MeV. 

the strength gained by the low-momentum modes. While it would be interesting to see if 
the inclusion of mode mixing in the quantal mean-field formulation reduces the enhance- 
ment, it will be hard to give a precise quantitative estimate of their effect. Therefore, 
it is probably more reasonable to consider the quantitative difference between the two 
approximations as an inherent uncertainty which can only be resolved by experiment. We 
note, however, that the above difference clearly demonstrates that a precise measurement 
of the dilepton spectrum not only can reveal the presence of DCC configurations but can 
also provide information about their size as well as their strength. 

3 A schematic model 

With the mean-field and semi-classical treatments, the time-dependent fields are dynam- 
ically coupled to the quasi-particle modes and the soft and hard modes of the fields are 
treated on an equal footing. To understand the main features of dilepton production due 
to the soft modes, we shall now consider the dilepton production process in a schematic 
model in which the soft and hard modes are coupled only via the electromagnetic interac- 
tion causing dilepton production. The hard modes are then represented by a thermal gas 
of pions having a specified temperature T, while the soft modes will be evolved dynam- 
ically according to the one-dimensional Bjorken expansion scenario. Hence the density 
operator can be factorized, p — pth® Pdcc, where pth is the normalized density opera- 
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tor for the thermal gas and pocc is the normahzed density operator for the soft modes, 
referred to as the DCC field. 

The current-current correlator W^j_p{q) can then be decomposed according to how many 
thermal pions are involved in the dilepton production process. The first term represents 
the coherent emission of dileptons from DCC fields alone and involves no thermal pions, 

Wl^J{q) = ldVy{Ux)jliy))ncce-^'-^'-'^ = {Uq)PM)Dcc ■ (24) 

where the DCC expectation value is (■ ■ ■)dcc = T^t^Ipdcc " " "]• 

Contributions to the current-current correlator from processes involving a thermal 
pion are 

Wl,]){x,y) = {7f^{x)n{y))Dcc{d^,n{x)d^7r\y))th- {d^rc^{x)d^7f{y)) Dcc{T^{x)rf\y))th 

+ {df,Ti\x)Ti{y))Dcc(j^{x)dy'K^{y))th - {df,Ti{x)Ti\y)) Dcc{T^\x)dyTi{y))th 
+ c.c. , (25) 

where we have omitted the vanishing disconnected part and all other terms containing 
{'ii\x)Ti\y)) th or {Ti{x)Ti{y))th which vanish in a thermal equilibrium environment. Using 
Eq. (0)? "we obtain 

+ (l + n^)(7r(g + A;)7rt(g + A;)) 

■ [n+(7r(g - k)T,\q - k))Dcc + % (7rt(A; - g)7r(A; - q))DCc]] , (26) 

The first two terms in the above equation correspond to the emission of one pion together 
with a dilepton by the DCC field. The factor 1 + is a result of the Bose enhancement 
in the final state. The second two terms proportional to represent the annihilation or 
absorption of one thermal pion by the DCC field. 

From Eq. ( p6[ ) one can readily see how the momentum distribution of the DCC field 
is imprinted onto the dilepton spectrum. Once the momentum of the dilepton is large 
compared to the inverse of the DCC domain size, the integral no longer has support 
from the Fourier transform of the pion field {n^k ± q)%'{k ± q))Dcc, which restricts the 
contribution to low momenta and to a small window in invariant mass. For example, 
if one considers a classical field that oscillates with a typical frequency of u; ~ and 
has a Gaussian distribution of width i?^ in coordinate space, one can see from Eq. (|26|) 
that the dilepton yield will be concentrated around an invariant mass of M ~ 2171^^. The 
width of this distribution will be of the order of 1/R_l in invariant mass as well as in the 
three-momentum q. 

Within this schematic model, it is possible to take approximate account of expansion 
by subjecting the DCC field to a boost-invariant Bjorken expansion. The DCC field 
depends then only on the proper time r = \Jt'^ — and, in the linear sigma model, its 
equation of motion becomes 

+ -0 + A(0^-w^)0 = ifa. (27) 
r 
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The initial values of and are sampled from normal distributions with suitable width 
parameters. The corresponding energy-momentum tensor is then relatively simple, 



T'^" = exi^u^'u'' - ^'^") + ey^^^ (28) 

where the four- velocity is = x^/t with = {t,0,0,z). The contribution to the 
energy density from the time dependence of the field is ex = |0^, while the contribution 
from the interaction is ey = 1(0^ — f^)^ — Ha — egs, with the vacuum energy being 
Cgs = jifn ~ '^^)^ ~ -^/tt- The expansion causes the energy density to drop steadily, 



§^{e, + ey) = , (29) 

and a simple power behavior is quickly approached, e/^ + ey cqTo / r. The coefficient cq 
then provides a convenient means of characterizing the particular solution. The limiting 
1/r behavior of the energy density is a general characteristic of the boost invariant scenario 
IP and we therefore assume that the energy density of the thermal gas drops in the 
same manner. This is accomplished by using a time-dependent temperature for the gas, 
r=(ro/r)V4To. 

We can take approximate account of the finite transverse size of a DCC domain by 
giving the chiral fields a common profile (such that the orientation of the pion field in the 
Cartesian space only depends on the proper time), 

g{x^) = exp(-xi/2i?l) . (30) 

This transverse profile is simply an ansatz and does not follow from nor is subjected to the 
equation of motion. Such a simple Gaussian form was suggested by numerical simulations 
and it suffices for our present purpose. The electromagnetic current density then 
becomes [[T^ 

j,{x) = V, fl^ g\x^) = vs fl"^ exp(-xi/i?i) . (31) 

where the dimensionless normalization coefficient is determined by the specific initial 
condition. 

The longitudinal expansion causes the current density to decrease as the proper time 
grows. As a consequence, the contribution to the dilepton yield from the coherent dilepton 
emission remains finite. Following Ref. and averaging over the ensemble of initial 
conditions, we find 



dygdqldM 24 ^ 6 m^MI 



where the transverse mass is M_l = yM^ + q\_ (with = q^) and Jq and A^o are Bessel 
functions. This differential yield is independent of the dilepton rapidity because of the 
boost invariance of the DCC field, and so it depends only on the initial value of the current, 
as determined by the ensemble average -< vl >~. It may also be noted that the coherent 
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emission vanishes as the transverse momentum of the dilepton photon approaches zero. 
This further diminishes the importance of coherent emission compared to the incoherent 
emission at low transverse momenta where the effect of the DCC is expected to be largest. 

In order to calculate the incoherent dilepton production by the DCC field, we need to 
perform the Fourier transformation of pion field tt{x) = T^{T)g{x±). We find 



rdr 7r,(T) [-iVo(M^r) T ^MM^r)] , po > \p\\\ 
rdr 'Ki{T) Ko{M±t) , \Po\ < \p\\\ 



(33) 



where = ^\pl and (^(g^) = 7ri?5_ exp(-g^i?^/4). 

One can then use (7r'''(g)7r(g)) = |[(|vri(g)p) + (|vr2(g)p)] in Eq. ( ^61) to calculate the 
dilepton yield from the interaction of the thermal pions with the DCC field, 



dN, 



(1) 



1 



dyqdqj_dM 



M2 



+ 



Stt^ M J 2uk{2TTy 

{l + nt){7r\-q-k)7r{-q~k))DCC 
l + n,){n{q + k)n^{q + k)) 

DCC 



+ nt{Ti{q-k)Ti\q-k))DCC + nl(ji\k-q)Ti{k-q))occ]- (34) 

Due to the boost invariance, we need only consider dileptons with zero rapidity so that 
q = {M±, qx, 0). Given the simple Gaussian transverse profile of the pion field, one can 
carry out the azimuthal integral, and the remaining two-dimensional integral can then be 
evaluated numerically by use of the Fourier transforms in Eq. (^3]). 

In the numerical evaluation of the dilepton production from the DCC field, we have 
employed the following rather conventional parameter values, A = 19.97, v = 87.4 MeV, 

= 92.5 MeV, = 135 MeV, with which H = f^ml = (119 MeV)^ and = 

600 MeV. We first calculate the Fourier transformation of the pion fields 
3f) for the solution of the equation of motion. We then calculate 



2A/2 + ml 
according to Eq. 

(7r"''(g)7r(g)) averaging over the initial configurations which then is used for the numerical 
evaluation of the dilepton yield in Eq. (0). The ensemble average -< 
and used in the calculation of the coherent emission yield, Eq. (|3^). 

Figure |^ shows the differential dilepton yield from a total space-time volume 



I ;^ is also extracted 



/'"max ~ ''"o ) due to incoherent emission from the DCC field together with the corre- 



spending average thermal rate. We have used the following parameter values, To=145 MeV, 
Rx=2 fm, tq=1 fm, rinax=102.4 fm, and the characteristic energy density carried by the 
DCC field was eo=58 MeV/fm^, which is the same as the thermal pion energy density at 
the employed temperature Tq. We have also separated the incoherent dilepton rate into 
bremsstrahlung and annihilation (absorption) part, since they are caused by very different 
processes and should have different features in the invariant mass spectrum. 

We find that the coherent emission rate in Eq. (^) is negligible as compared to the 
incoherent emission. This is related to the fact that the coherent emission in our scenario 
depends only on the third component of the conserved isovector current. Because the 
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Figure 7: Dilepton spectrum from the DCC field, botfi bremsstralilung (dotted), annihi- 
lation and absorption (dot-dashed), and their sum (solid), as well as the thermal emission 
(dashed). The initial temperature of the thermal environment is To=145 MeV, while the 
characteristic energy density carried by the DCC field is eo=58 MeV/fm^. 



energy density carried by the isovector current decreases much faster than that carried 
by other degrees of freedom, the associated dilepton emission is thus also less important. 
However, the incoherent dilepton emission is very significant in comparison with the ther- 
mal production. The spectrum from the pion annihilation and absorption by the DCC 
field (dot-dashed) has a structure manifest of two components in the second two terms in 
Eq. ([3^), depending on the energy flow. In the contribution from the annihilation with 
the DCC field, energy flows out of the DCC field, therefore the dilepton spectrum has a 
threshold at M = 2171^^, if the DCC pion field oscillate with a minimum frequency of m^r. 
If the energy flows into the DCC field, dileptons are then emitted when thermal pions are 
absorbed by the DCC field. The dilepton spectrum from these absorption processes has 
no threshold and dominates the annihilation and absorption spectrum at small invariant 
masses M < m^, as seen in Fig. |^. 

Shown as the dotted line in Fig. |^ is the dilepton spectrum due to pion bremsstrahlung 
from the DCC field, the first two terms in Eq. (0). Since the emitted pions have at least 
a minimum energy of m^, the DCC field must have some higher frequencies in order to 
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emit a thermal pion plus a pair of leptons. This means that the DCC must have some 
quasi-particle modes whose masses are larger than the pion mass m^. To illustrate this, we 
show in Fig. ^ both the pion and sigma fields as functions of time and their corresponding 
Fourier spectra. As one can clearly see, besides the normal pion mode with mass m^^ (this 
is the frequency of the fields for the zero-momentum mode we are considering), there are 
others resonances with higher masses. While the resonances at ~ 3.5m^ and ~ 5.5m^ are 
clearly visible, the higher ones cannot be seen here due to their small amplitudes and the 
limited resolution due to the finite time interval used in our numerical calculation. These 
quasi-particle modes are normally referred to as parametric resonances in a non-linear 
and strongly coupled system. We will not elaborate on the interesting physics associated 
with the parametric resonances, but just point out that they play an important role in 
the preheating of the early universe and the amplification of the long-wavelength 



pion mode following an extremely nonequilibrium initial condition |T8|, ^ With these 
parametric resonances in mind, we can readily understand the spectrum of the dilepton 
from the bremsstrahlung processes: dileptons are produced via the transitions of the high 
resonances to the normal pion mode. As we increase the invariant mass of the dilepton, 
some of the transitions are gradually turned off as the energy of the dilepton becomes 
larger than the mass differences between the higher resonances and the normal pion mass 
rriTt. This is why there is a slight oscillation in the bremsstrahlung dilepton spectrum. The 
oscillation vanishes when the invariant mass reaches the highest parametric resonance that 
the numerical solution can resolve. This is also the case when we increase the transverse 
momentum of the dilepton. 

In Fig. 1^ is also shown the sum of the different contributions to dilepton spectrum 
from the DCC field (solid curve) as well as the contribution from annihilation in the 
cooling pion gas. In general, one can see that the incoherent dilepton production below 
and near 2m.,^ threshold region is significantly larger than the thermal production. As 
we already pointed out, because of the finite spatial size of DCC field in the transverse 
direction, dileptons from the DCC field exhibit a much faster decrease with the transverse 
momentum q± than those due to thermal production. 



4 Summary 

We have calculated the production of dileptons from disoriented chiral condensates using 
a quantal mean-field as well as a semi-classical treatment for the time evolution in the 
linear sigma model. We have compared the dilepton spectra obtained when using so called 
quench initial conditions, which lead to a strong enhancement of long wave length pion 
modes (DCC), with those obtained from thermal initial conditions. Compared to the 
thermal spectrum the quench initial conditions lead to a strong enhancement (factor 20 - 
100 depending on the model) at an invariant mass of about M ~ 2m,r- This enhancement 
is confined to dilepton momenta of g < 300 — 500 MeV and also rather narrow in invariant 
mass. Based on the analysis in a schematic model we find that the observed enhancement 
at M ^ 2m7r is due to the annihilation of thermal pions with those from the DCC. 



18 




-0.1 







25 50 75 100 

T (fm/c) 

0.15 0.15 




> 

CD 

CD 



25 50 75 100 

T (fm/c) 



2 4 6 8 

co/m„ 



Figure 8: The time evolution of the tti and a DCC field and their corresponding Fourier 
spectra. 



We furthermore found some evidence for an enhancement at low invariant masses, 
M < mjr, which, however, we could not unambiguously establish due to numerical back- 
grounds. In the schematic model this enhancement is due to bremsstrahlung-type pro- 
cesses involving the absorption or emission of pions by and from the DCC. The latter 
reflects the rich dynamical structure of the DCC field. The spectral distribution of the 
pion field shows distinct peaks not only close to the pion mass but also at higher frequen- 
cies. This phenomenon is related to the so-called parametric resonances and is responsible 
for a comparatively strong dilepton yield as a result of pion emission processes. 

Within the schematic model, we could also address the more realistic scenario of a 
longitudinal expansion. We find that the above enhancement remains also in this case. 

Both dynamical solutions, which are based on distinct approximations to the linear 
sigma model, predict qualitatively the same enhancement: A large bump at M ~ 2771^^ as 
well as an enhancement at low invariant masses. Overall, the quantal mean-field model 
seems to yield a larger enhancement than the semi-classical treatment. This quantitative 
difference may be due to several factors. First, the Bose enhancement factors are absent in 
the semi-classical treatment. Second, the semi-classical treatment incorporates the mode 
mixing resulting from the non-linear interaction and this mechanism tends to reduce the 
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number of pions in the DCC state. Both of these features lead to a somewhat smaller 
signal. 

As far as experimental observation of this enhancement is concerned, it will be prob- 
ably very difficult to see the enhancement at low invariant masses, M < m^, because 
this region will be dominated by the Dalitz decay of the n^. The enhancement around 
M ~ 2m^ on the other hand should be observable in principle. In this mass range the ma- 
jor competing channel is the Dalitz decay of the rj. From the analysis of CERN SPS data 
(see e.g. [^) the 1] Dalitz is at most a factor of five stronger than the pion-annihilation 



channel. Therefore an enhancement of the pion annihilation by a factor of ten or larger 
should be visible. However, due to the background from false pairs it is very difficult to 
measure dileptons of small momentum where the expected enhancement is located. 

In the present work we have ignored expansion, with the exception of the schematic 
model. This will have to be taken into account in future work for a more realistic de- 
scription of a heavy-ion collision. Morevoer, expansion provides one possible justification 
of the quench initial conditions [0. Certainly in a given effective theory, a consistent 
treatment should generate the quench conditions dynamically starting from a more or 
less thermal state at temperatures above Tc. The resulting dilepton spectrum from such 
a calculation should then provide a more realistic picture of possible enhancements to be 
expected in heavy ion experiments. In addition such a calculation needs to account for 
additional dissipative processes which are not included in the linear sigma model, such as 
the scattering with vector mesons. 

However, too little is known about the physics that reigns close to the chiral phase 
transition to allow us to give a final theoretical answer about the existence of DCC states 
in relativistic heavy ion collisions. Therefore, experiments measuring unique observables, 
such as the one addressed here, will be needed to determine that issue and its implication 
for the chiral phase transition in matter. 
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